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ABSTRACT 


After surveying electrically-driven loads used on productive farms# the 
investigators chose three pumps for testing at voltages and frequencies far 
outside the normal operating range. These loads extract and circulate water 
and move heat via air# and all are critical to farm productivity. The object 
was to determine the envelope of supply voltage and frequency over which these 
loads would operate stably for time intervals under 1 hour. This information 
is among that needed to determine the feasibility of supplying critical loads# 
in case of a utility outage# from a wind-driven alternator whose output vol- 
tage and frequency will vary dramatically in most continental wind regimes. 

The report surveys other related work, describes the salient features and 
limitations of the test configurations used# describes the data reduction# and 
outlines the development of simulation models suitable for a small computer. 
The results are primarily displayed on the voltage- frequency plane with the 
general conclusion that the particular pump models considered will operate 
over the range of 50 to 90 Hz and a voltage band which starts below rated# 
decreases as frequency decreases# and is limited on the high side by excessive 
motor heating. For example# centrifugal pump operating voltage ranges as 
extensive .4 to 1.4 pu appear possible. 

I&rticular problems with starting, stalling due to lack of motor torque# 
high speed cavitation# and likely overheating are addressed in a listing of 
required properties for wind-driven alternators and their controllers needed 
for use in the isolated or stand-alone configuration considered. 
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I. INTRODUCTION AND PROBLEM STATEMENT 


Before reliable utility electric service was widely available on farms, 
non-vehicular mechanical power was provided by other sources including mecha- 
nical hydro power for grinding, wind-driven pumps for water supply, and draft 
animals harnessed to lift water for irrigation. Reliable electric motors 
displaced these sources and have made a significant contribution to agricu- 
ltural productivity. However, there is a small possibility of long-term 
utility outages, due to fuel shortages or disaster. Could alternative sources 
then continue farm production when utility electric energy is unavailable or 
too expensive to deliver to remote sites? This report is a preliminary step 
towards the answer. 

Although a complete list of electrical farm equipment is extensive, there 
are several motor loads necessary to almost all farming operations. They 
include water pumping and circulation for animal, crop, and residential needs. 
Cooling and heating, while not necessary in some farm operations, is needed 
for rural residences as well as for dairy operations. 

The alternative source considered in this study is wind energy as conver- 
ted by variable-speed self-excited alternators — a practical wind energy alter- 
native not requiring utility power for excitation or control. The common 
wind-driven induction generator is currently practical only when excited by a 
strong outside source. Photovoltaic sources would require a force-commutated 
DC to AC inverter designed for isolated operation since farm motors are all AC 
and designed for nominal utility service voltages and 60 Hz. Fueled auxiliary 
generators can provide near utility-grade energy and present few new technical 
problems for evaluation when the load is nearly constant.. In addition, 
during a prolonged utility outage, that fuel could instead be used for farm 
vehicles if wind energy was available. 

A stand-alone wind-driven alternator, without extensive and expensive 
controls, will vary unpredictably and substantially (4 to 1) in speed (and 
thus frequency) and in terminal voltage as wind velocity and load demand vary. 
Thus, any load connected to such a source must operate over a voltage-fre- 
quency range (envelope) far in excess of normal load design conditions. The 
latter are usually 60 Hz + 5%, and 120, 240, or 480 volts + 5%. Although AC 
motors can be manufactured to operate over a wider range of input conditions, 
the cost may be excessive given the low risk of an extended utility outage. 

In other words, why retrofit existing electric drives with expensive wide- 
range motors for use with a new windmill when the extra range is needed only 
for a low-probability prolonged outage? Perhaps it is better to augment the 
controls of a wind-driven alternator so its output can drive standard motors 
designed for small variations in frequency and voltage under the much greater 
voltage and frequency deviations present in isolated operation. 

The essence of this problem is: Find the operating envelope in voltage 
and frequency of several critical motor-driven farm loads and find the proper- 
ties of the controller of an alternator windmill to enable it to operate these 
loads outside design conditions without immediate damage. Although a motor 
and load may appear to operate satisfactorily for a short time outside the 
design envelope, extended operation under such conditions may reduce service 
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life and increase maintenance costs. For example, temperature rise and vibra- 
tion limits may be exceeded. If extensive off-design operation is antici- 
pated, these conditions should be investigated. A thorough exploration of 
this problem was far beyond the time and resources available under this grant 
so the effort was simplified. 

A useful first step is finding the short-term (several minute) steady- 
state operating envelope of three common and commercially available farm 
pumping loads for frequencies between 40 and 100 hz and for voltages 
ranging across and far outside the manufacturers specifications. 

A field feasibility study, including limited testing, as opposed to a control- 
led lab investigation, was intended. There were some limitations. For ex- 
ample, shaft torque and speed were impossible to measure directly due to pump- 
motor packaging and location charateristic of these commonly used pumps. 
Computer simulations were then used to compare with test results and together 
they suggest windmill controller characteristics for isolated steady-state 
operation of these critical loads. Future work is needed to investigate 
service life, transient conditons like motor starting and sudden load changes, 
and the effect of these loads on the behavior of wind turbine generators. 
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II. RELATED WORK 


Many modern textbooks on electric machines [eg: Ref. 1 and 11] include 

development of circuit and dynamic models applicable to caged induction 
motors operating over a wide range of supply voltage and frequency. The 
emergence of a substantial market for solid-state variable-speed induction 
motor drives is testimony that motors can successfully operate over a wider 
range of supply conditions than is suggested by the motor nameplates. The 
specification sheets provided by manufacturers for various classes of caged- 
rotor motors list the effects of 10% voltage and 5% frequency changes on motor 
torque r current, speed, power factor, efficiency and temperature rise. For 
example, one line of integral horsepower drip>-proof motors is said to decrease 
torque by 10%, increase speed by 5%, and "slightly increase" efficiency for a 
5% frequency increase. 

Although a literature search was not part of this study, the investiga- 
tors are not aware of an extensive literature reporting on large changes in 
supply parameters for farm pumps except for a 1983 USDA study [Refs. 1,12]. 
This involved laboratory testing with a permanent magnet three-phase 9 kw 
alternator designed for a windmill feeding a centrifugal pump. The instrumen- 
tation used was more complete and accurate than that available for the field 
tests reported here and the permanent magnet alternator restricted motor 
voltage control. Nevertheless, the investigators found that "standard three- 
phase, AC induction motors will pump water when operating at 30 hz and 85 
volts" when the motors were rated at 230 volts and 60 hz. Their work charac- 
terized motor power output, efficiency and temperature rise in a complete 
manner not possible in the tests reported here due to inaccessibility of the 
motor shafts in commercial pumps. 

Direct water pumping and heating is also possible using wind energy. The 
use of wind-driven centrifugal pumps for water heating has been explored [Ref. 
3] and this method of heating could be compared with a heat pump driven by a 
stand-alone electric windmill as reported in this study. Models of pump 
behavior are widely available [Ref. 4] and the behavior of wind-driven pumps 
is described in a recent book [Ref. 5] . 

The control of a separately-excited wind driven alternator has been 
addressed by a number of manufacturers of windmills. Some recent innovative 
work is described in [Ref. 13] . 

The behavior of single-phase submersible pumps was not referenced in 
these papers and neither was work on single- phase heat pumps. Nevertheless, 
the insights presented in the referenced papers aided the investigators in 
this study and verified several of the required controller properties listed 
in Section VIII. 
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III. CHOICE OF TEST LOADS 


A generic listing of electric farm equipment is shown in Tbble 1. Most 
devices are rated at 115 or 230 volt single-phase 60 Hz since three phase 
service is not available on most farms or if available, is more expensive than 
single-phase service. Multiple single-phase farm services totaling 800 amps 
at 240 volts are said to exist. The National Rural Electric Cooperative 
Association (1800 Massachusetts Avenue, Washington, DC) is one source of 
agricultural load information and the US DA at Ames, IA is another. A study 
sponsored by NRECA [Ref. 6] reports peak demand for a 120 cow dairy farm as 29 
KW with about 1/2 of this load due to electric motors. The brief survey done 
for this report found most motors at 230 volts, single phase with nameplate 
horsepower mostly under 7.5 hp but with some motors purportedly as high as 50 
hpt The majority of motors are fractional or 1 to 5 hp motors with caged 
rotors. 

Many materials handling applications have high starting loads and are 
used intermittently such as conveyers and loading and unloading of biomass 
materials. 

Ventilation applications (fractional to 2 hp) often pump a fixed air-flow 
volumetric rate against a small positive pressure for air supply and tempera- 
ture control of animal housing (hogs, poultry). Grain drying and aeration 
operations require sufficient pressure to penetrate piled grains. Some ma- 
terials, such as chopped silage, is moved by air flow so velocity in such 
systems must be maintained to silo heights. 

Water extraction, circulation and heating is done on every farm and since 
storage is possible in ponds and tanks, can be done intermittently when alter- 
native power is available. 

Refrigeration is necessary in dairy operations and poses special problems 
due to high starting loads and to the sealed motor-compressor configurations 
usually used. Heating is used for crop drying and animal housing in cold 
climates. 

The investigators chose for testing and analysis three of these loads 
which they believe are critical to many agricultural operations. These are: 

* The extraction of water from deep wells. 

* The circulation of water from storage reservoirs for animal watering 
and crop irrigation. 

* The heating of air by a heat pump since it involves both air 
circulation across the condensor and evaporator and a hermetic 
motor-compressor package as used in refrigeration. Attempts were 
made to borrow a milk cooler to represent a hermetic cooling load 
but only an air-to-air heat pump-air conditioner could be obtained 
and tested. 
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FEED PREPARATION — 


WATER SYSTEMS — 


buir milts 

- - - - • - hjmmtr mills 

roller mills 
roughage shredder 
feed blenders 
feed meferer 
feed mills, automatic 
feed mill equipment and accessories 
corn shelters 
seed and grain cleaners 

FEED HANDLING AND CONVEYING EQUIP. — 

farm scales 
silo unloading devices 
bunk feeders 
feed conveying equipment 
bucket type elevators 
hay mow conveyors 
sweep un-loader 
grain bin un-loader devices 
wagon and truck un-loading devices 
wagon loading device 
auger or scrow conveyors 
elevators and/or conveyors 
fruit and vegetable equipment 

CROP AND GRAIN ORYERS — 

batch dryers 
continuous type dryers 
heated air units 
supplemental heat units 
aeration fans 
crop moisture meters 
crop drying accessories, bins 
grain probe thermometer 


ENVIRONMENTAL CONTROLS — 

fans, ventilating 
walk-in coolers 
insect destroyers 
soil pasteuriser 
insulating glass 
switches 
thermostats 


EMERGENCY GENERATING EQUIPMENT — 

electric generating sets (gas) 
electric generating sett (diesel) 
electric generating sets (tractor driven) 
motors and accessories 

ELECTRIC HEATING — 

electric heating 
heating cable 
temperature controls 
in-line heaters 
crop drying equipment 
environmental control 


Table 1 Farm Electrical Loads 


pumps (shallow well) 
pumps (deep well) 
pumps (portable) 
hot water heaters 
water softener systems 
wafer purification 
water systems accessories 
irrigation system equipment 

MANURE DISPOSAL — 

barn eloanort or gutter cleaners 
liquid manure disposal systems 

DAIRY EQUIPMENT — 

milking machine units 
pipeline milking units 
bulk coolers 
walk-in coolers 
chums 
pasteurizer 
heat recovery unit 
dairy rontilation 
heating units 
hot water heaters 
dehorners and branders 
clippers 
dairy housing 
cattle cleaner 


POULTRY EQUIPMENT — 

incubators 
brooders 
poultry feeders 
poultry waterers 
egg cleaning and handling equipment 
laying nests 
poultry ventilating equipment 
walk-in coolers — other coolers 
fountain medicstor 
debeaker 
sealder 
time controls 
poultry housing 

FEEDLOT EQUIPMENT — 
CATTLE, SHEEP, SWINE 
EQUIPMENT AND HOUSING 

self feeders 
waterers 
ventilation equipment 
water level control unit 
heating units (floor heaters) 
heating units (tank daicer) 
de-horners and branders 
clippers 

electric fence controllers 
gates 
farm scales 

SHOP EQUIPMENT — 

air compressors 
welders 
soldering tools 
power tools 
cement misers 
battery chargers 
steam cleaners 
heating units 
engine starting units 
incinerators 
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The additional requirements included single- phase 60 Hz units in common 
use (mass-produced) with horsepower under 2 [ 1.5 kw] so that the 5 KVA power 
source used for the tests could start them under load, and permanent split 
capacitor windings so there were no centrifugal- switches to affect low speed 
and thus low frequency operation. This is a significant restriction since 
many motors use internal centrifugal starting switches to disconnect a star- 
ting winding at about 75% of normal speed. 

Table 2 summarizes the specifications of the tested loads. The investi- 
gators thank the manufacturers of these devices for providing the further 
design information necessary for computer simulation of the behavior of their 
products. 

Since the project was planned to evaluate feasibility of off-nominal 
operation in agricultural applications, the devices were field tested under 
conditions that were intended to represent typical applications. While the 
feasibility of operating other farm loads with wind power could be tested, the 
investigators believe that water pumping operations are necessary for all 
farming operations. Heat pumping to and from air is less critical for non- 
dairy operations but is of importance for space conditioning for residences as 
well as for farm building and so was included. 

The submersible pump was suspended in a steel well casing inside a well- 
house with the plastic output pipe feeding a well-head flow meter and pressure 
gage before discharging thru a gate valve into a long plastic pipe. This 
setup is diagrammed in Fig. 1 and a photo shown in Fig. 2. 

The circulating pump, a plastic self-priming rotor inside a cast housing, 
drew water from a farm pond thru a foot valve and short pipe and its output 
flow fed thru the same flow-meter, pressure gage and valve as above before 
entering a plastic pipe 110 meters long which returned the flow to the pond at 
a small head. The set up is illustrated in Fig. 2 and 3. The initial test 
used a stock tank and a shorter pipe before the pond was used. 

Hie air conditioner-heat pump was a typical 1 1/2 ton residential unit 
with an exterior enclosure containing controls, hermetic motor-compressor (the 
latter with plug pistons), shaded-pole fen and outside heat exchanger. Two 
sta n dard-lengths of precharged refrigerant lines fed the inside heat exchanger 
(baffled "A" coil). This coil divided a two-chamber plywood "furnace" with a 
squirrel-cage fen mounted at the bottom. The fan forced air upward thru the 
"A" coil, the upper plenum and out thru an opening 30% larger than the fan 
inlet area. The "furnace" (see Fig. 4) was instrumented with thermistors and 
draft gages. Hie seb-up is illustrated in Figs. 5 and 6. This unit was 
loaned to the investigators by its manufacturer. Without this loan, the 
project could not have been completed. 
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Property 

Submersible 

Pump 

Circulation 

EUmp 

Heat Outdoor Furnace 

Pump Fan Fan 

Motor Motor 

Variable 

Supply 

Alternator 

Manufacturer 

F.E. Myers Co. 
Ashland, OH 
44805 

F.E. Myers Co. 
Ashland, OH 
44805 

Addison Products Inc. Dayton 

Addison, MI Elect. 

49220 Mfg. Co. 

Qiicago 
IL 6009! 

Miller Elect. 
Mfg. Co. 
Appleton, WI 
54912 

Model Number 

J10-18 

QP10 with 
B828 motor 

QH818 4C684 

AEAD-200LE 

Fated Freq. Hz 

60 

60 

60 60 

60 (100) 

Rated Voltage 
RMS 

230 

230 

230/208 115 

240 

fated Amperes 
RMS 

7.9 

mm 

9.3/10.0 8.9 

w/o fan Cl. 9) 

5 KVA 

Locked Rotor 
Amps 

34 

34 

48.4 4. 4 

44 (short 
circuit) 

Run-start 

Capacitance 

25 uf 

warn 

25 uf shaded 

pole . 

. 

Rated horse- 
power 

1 

i 

1/4 1/4 

16.5 of gas engine 
prime mover 

Rated speed 
rpxn 

3430 

3450 

2.075 1050 

1800 (3000) 

Rated Pump 
Efficiency 

51.8% e 
16 gpn 

63.84 8 
40.5 gpn 



Motor Service 
Factor 


1.4 (100% duty 
cycle 40° C) 


100% duty cycle 

Rated volumetric 
flow (near 
actual operating 
conditions) 

22.5 gpn 8 
115 feet 
total head 

19.2 gpn 

2670 1320 

cfm fan 


Rated Fewer 
Output 



19,300 BTq/hr 
8 47° beating 
17,800 BTQ/hr 
8 47°F cooling 


Rated Motor Input 

1.22 KW 


2.17 KW for above 


Load During 
Tests 

100 ft 1 inch 
I O pipe to 
instrument 
manifold plus 
150 ft outlet 
pipe 

20 ft. 1 1/4"D 
inlet pipe to 
manifold + 25 ft. 
1 1/4"D pipe to 
300 ft 1*D pipe 

1 

Ambient air 
59 - 71 °F 



Table 2 - Specifications and Properties of 
Tested Loads and Source 
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Schematic of Submersible Pump Test 



Figure 2 

Photos of Water Pump Tests 
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Schematic of Circulating Pwrp Test 











Figure 5 

Photo of Heat Pump Test 
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170-400 VAC from electrical 
instrment package 



Figure 6 


Schematic of Heat Pump Test 





IV. MEASUREMENT TECHNIQUES, LIMITATIONS AND DATA REDUCTION 


1) lOStCUHEOtS : 

A sophisticated computer-based instrumentation system for this project 
would have absorbed the entire budget. Because the problem was simplified to 
determination of the steady-state operating envelope including power input and 
output variables, it was possible to use classic moving- pointer instruments 
and inexpensive field-grade transducers ordinarily used for commercial appli- 
cations. The electrical input measurements were made with nearly identical 
setups for all three applications. Moving iron meters were used for voltage 
and current and a moving coil meter for power. Voltage and current readings 
were occasionally checked with a digital volt meter and clamp-on ammeters 
while frequency was measured with a digital meter and checked with an analog- 
output frequency to voltage hybrid integrated circuit package. A digital 
tachometer was used to verify the generator frequency measurements and to 
measure shaft speed on the furnace fan — the only load with an exposed shaft. 
Water flow was measured with a conventional commercial watermeter with a 
pulsed electrical output. Water pressure was read from a commercial Bourden- 
tube gage calibrated in psi. The flowmeter was checked by timing the filling 
of a 5 gallon can. The pressure measurements were not redundant — an unfortu- 
nate omission because pressure measurements are used to calculate pump power 
output and so should be double-checked.. For the heat pump test, temperatures 
were measured by thermistors with a digital readout and checked by a thermo- 
couple — which was used with a wet sponge for wetbulb readings. Relative 
humidity was measured with a cheap "hardware store" gage. Static pressure in 
the "furnace" was measured with an inclined draft gage and roughly checked 
with a water-filled U-shaped manometer tube. Air velocity was measured with a 
mechanical pickup tube and meter. Load current or voltage waveforms were 
observed via a wideband toroidal current transformer or a toroidal potential 
transformer into a battery-powered oscilloscope. An extra set of thermistors 
was used to record the transient behavior of the heat pump piping temperatures 
to estimate time constants in the refrigerant loopt No flow or pressure 
transducers were available for monitoring the refrigerant loop. Power for the 
instruments that needed it in the water pumping tests out in the field was 
provided by a vibrator inverter of indeterminate age but impeccable relia- 
bility. 

The test setups are diagrammed and photographed in Figs. 1 thru 8. The 
instruments and transducers are listed in Table 3. 


2) Test pCQtQCQl: 

After the test setup was checked, the gasoline- powered welding generator 
was started and allowed several minutes warmup time. This generator was 
chosen because it will operate over a wide speed, and thus frequency, range 
without damage to its windings and because its exciter circuit allowed for 
considerable output voltage variation. However, the voltage variation was not 
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Figure 7 
Instruments 
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Schematic of Electrical Measurements 
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Table 3 

Instruments Used 
in Tests 



as wide as the investigators desired and was increased for the heat pump tests by 
two 13 kva autotransformers operated in parallel. Furthermore, the generator 
governor was modified to allow steady^state operation over a range of speeds 
corresponding to 40 to 100 Hz. 

With the generator warmed and operating at 60 hz and 230 volts, the load 
was switched on thru a fused switch and allowed to stabilize.. This took 
about 15 minutes for the heat pump but only a few minutes for the water pumps. 
After several measurements at nominal conditions, the generator throttle, 
field control, or autotransformer was varied to obtain a new test condition. 

The test conditions obtained and the time are noted in the data sheets repro- 
duced in Appendix A from transcription of the rough field notes of the princi- 
pal investigator. At the edges of the operating envelope, operation was 
unstable and the load motors either stalled due to inadequate torque, the 
water pumps cavitated, or the motors took so much input power that the run was 
terminated after or during the taking of readings. Motor overheating was 
noted only once when the heat pump motor thermal cutout opened and stayed open 
until the hermetic container was cooled. No motor winding temperatures were 
observed due to inaccessibility of the windings. After any shutdown, several 
minutes passed before a restart was attempted. For the heat pump, a longer 
interval was needed for the refrigerant pressures to equalize. 

Approximately 170 sets of readings were taken during the successful tests 
over a period of several months and they all are included in Appendix A. The 
test site climatic conditions were moderate for the circulating and heat pump 
tests but below freezing for the submersible pump test which affected only the 
investigators and the instruments because the submersed pump was near 10 tf C. 


3) Test limitatiQos: 

Inspection of the data sheets and the previous paragraphs suggest several 
limitations of the data — which the investigators believe do not change the 
conclusions of the study. However, if the study was to be repeated or exten- 
ded, they should be considered in future test planning. 

For several heat pump runs at high frequency, the measured real power 
slightly exceeded the calculated apparent power — which is physically 
impossible. Wattmeter error is suspected since the current and 
voltmeters were rated from 25 to 500 Hz while the wattmeter is not. 

The digital readout for furnace input and output temperatures yields 
integral numbers so temperature differences across the "A" coil 
could be in error by as much as 7% with a corresponding error in 
heat power output. 

The flowmeter measurements were timed with a watch and small errors 
are likely. 

The load voltage was measured at the instrument package, not at the 
motor terminals at the end of 3 to 30 meters of #8 or #12 wire. 

This error might be as much as 2%. 
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* The generator governor could not always hold the frequency steady so 
the frequency regulation had a total variation (peak to peak) of one 
Hz. 

' The pressure drop across the flow meter was estimated from the 

manufacturers data and not measured. Similarly, the pressure drop 
across the foot valves and the gate valve and the pipe couplings 
were estimated — not measured. This causes errors in the load 
characterization. 

* The pressure gage had no negative pressure calibration marks so the 
circulating pump suction was estimated. Hus is likely to cause 
pump power output errors as much as 10%. 

In spite of these errors, the investigators found the data set in general to 
be consistent and relatively close to values predicted by the computer simula- 
tion. 


4) Test data reduction: 

Calculations using test data were done for two reasons. The first was to 
describe the pump electrical input behavior in terms of input kva (apparent 
power) and power factor as well as pump output in kW. Hie latter required 
complex calculations and corrections since the mechanical outputs of the 
submersible and heat pumps were not directly available for measurement. The 
other reason was to model the pump loads so they could be BBincluded in the 
computer simulation models to allow comparison between predicted and measured 
performance. Hie results of the data calculations are entered in the columns 
marked by an asterisk (*) in the data sheets which are in Appendix A. The 
investigators regret the use of English units but their instruments were so 
calibrated and excessive conversions increase the possibility of error. How- 
ever, the input and final output variables are presented in SI units. 

The data reduction for the electrical input variables required few assump- 
tions. The magnitude of the apparent power is just the pump input amperes 
times the pump input voltage. Hie burden of the electrical instruments was 5 
VA or less so no loading corrections were necessary. Hie ratio of measured 
real power from the wattmeter to the calculated apparent power gave the power 
factor of the pump under test for the input conditions observed. 

The pump power output calculations used data from less accurate gages and 
required more assumptions and are thus less accurate. However, since the same 
assumptions and resultant corrections were used for all runs of each of the 
three loads, the pump behavior is fairly characterized by comparing the power 
output behavior at nominal or rated conditions with the output under signifi- 
cant deviations from those rated or nominal conditions. 

Ecwet output corrections foe submersible pump: 

Since power is force times velocity or pressure times volumetric flow, 
and since the instruments were calibrated in English units, the S.I. formula 
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N 


/ nr^ 


Pcwer= Pressure --- x Flow — watts 


m? 


\sec/ 


is converted to 


Power = 0.4350 x GPM x PSI watts 

where GPM is flow in gallons per minute and PSI is total head in pounds per 
square inch — which are the units of the instruments and of the pump specifica- 
tions. If pressure had been measured at the outlet pipe discharge, and con- 
verted to head, then total head in this equation would have been this pressure 
head plus the difference in altitude between input and output water levels. 
This was not done so the following corrections were used to estimate pump 
mechanical output from the measured data. In the submersible pump measurement 
configuration used, this formula yields power measured at the surface manifold 
where the pressure is measured after the flow transducer. Corrections are 
needed for the pressure drops caused by: the flow transducer; the submerged 

foot or inlet valve at the pump input; and the friction loss of the 100 feet 
of pipe between the submerged pump outlet and the manifold. The power needed 
to lift the water flow against gravity is also required to estimate the pump 
output at its outlet. 

The power required to lift water against the gravity field is 
proportional to the net head times the volumetric flow. Since the net head 
(the difference in elevation between the water level in the well casing and 
the manifold at the surface) is about 7.93 m, the power required is 

P-Qf t = 4.904 x GPM watts 

Although derived here from first principles, this formula yields the same 
results as Johnson [Ref. 5, p. 282]. The pressure drop correction due to 
friction from the 100 feet (30.48 m) of 1 inch ID (.0254 m) smooth plastic 
pipe first requires an estimate of the Reynolds number (Rg) of the pumped 
fluid. Across the range of measured flows (about 5 to 30 GPM), R^ ranges 
from 1.38 x 10 4 to 8.27 x 10 4 . 

The relation is, for water at 50°F (10°C) , 

Rg = 2758 x GPM 

Assuming that the plastic pipe is very smooth, the corresponding friction 
factors range for f(R e ) of 0.027 to 0.018 [Ref. 7, p. 95]. The head loss 
Ah due to this friction is 



where z/D is the length to diameter ratio of the pipe, g is the gravity 
constant, and V is the flow velocity. In the units of the instruments, this 
becomes for this test configuration 
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Ah = f(GEM) x 3.109 X (GEM) 2 feet 

Since the pressure loss and the head are related by the fluid density, the 
correction becomes 

APSI = 1.347 x f(GPM) x (GPM)^ pounds per square inch. 

Thus, for each flow measurement, the corresponding pressure drop due to 
friction can be found. 

The pressure drop for the flow meter is given in the manufacturers data 
sheet and is included in Table 4. The pressure drop in the pump inlet 
valve is assumed equal to that of a half-opened gate valve of the same 
diameter or with K = 5.6 in the formula 

V2 


In terms of pressure loss, this becomes 

APSIfQQj. va ^ ve - 0.006288 (GEM) ^ pounds per square inch 

These pressure drops were added, converted to watts, and all corrections 
were added to give a total correction to measured surface output in watts as 
shown in Thble 4. The corrections are plotted in Fig. 9 which was used to 
estimate the correction for each pump run and added to the measured surface 
power to get the estimate of pump power output shown in the data sheets and 
used to illustrate the measured behavior of the pump in Section V. 


Eqv ec output fee© circulating pu©p: 

There was no correction made for the output calculation of the 
circulating pump. The output formula used was 

Pcwer = 0.4350 x GEM x ESI watts 

as in the case for surface power for the submersible pump. The suction 
pressure calibration was poor so these readings likely induce more error for 
this load where the PSI reading should be the pressure rise across the pump 
since the input and output water elevation was the same. 


Output calculations fee tbs beat pump: 

The thermal output power calculation from the furnace output vent 
driven by the outside compressor via the refrigerant loop and heat exchanger 
(the "A" coil) in the furnace requires several variables. They are the air 
mass flow (AMF) , the temperature difference AT across the "A" coil, and the 
specific heat of the air flow thru the furnace. Again, the instruments and 
heat pump components used were calibrated or specified in English units so 
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TSible 4 Submersible Pump Output Power Corrections 




these were used for the intermediate calculations although the output was 
expressed in kilowatts. 

The volumetric air flow is needed to find the air mass flow and the 
former was found from the furnace blower fan specifications and roughly 
checked by an air velocity profile taken at the exit vent to the furnace 
output plenum chamber. . The fan specifications (the Dayton fan is described 
in feble 2) give a flow of 1320 ftvmin against a head of 0.6 inches of 
water with a fan motor current of 8.9 amperes at 115 volts 60 Hz AC and a 
fan shaft speed of 1050 rpm. On the first day of testing, the measured 
static pressure in the fan plenum was 0.61 inches of water at 1047 rpm with 
a measured draw of 8.9 amps at 117 volts and 60 Hz and with a measured power 
input of just under 600 watts. Interpolation of the fen catalog data 
suggests 1313 ftvmin of flow. A few small leaks in the plywood furnace 
were taped before the second day of testing on May 21, 1985 and the pressure 
increased to 0.62 inches of water with 9 amps of current. The pressure 
increase suggests by interpolation a flow decrease to 1306 ftvmin but the 
increased current and a slight increase in fan power input suggests more 
flow so 1313 was used. A crude check was made by measuring the air velocity 
profile at the furnace exhaust vent with a velocity gage. The measured 
profile peaked at 2200 ft/min at the exhaust vent center and decreased to 
1900 to 2000 ft/min at a distance of 1/2 inch in from the edges of the vent. 
The flow velocity was fairly regular due to the diffusing effect of the "A" 
coil fins. The assumption of an average outlet velocity of 2000 ft/min over 
an outlet area of 8 by 12 inches yields a volumetric flew of 

8 x 12 

---- — x 2000 ft/min = 1333 ftvmin, 


which checks the 1313 figure interpolated from the iranufacturers data. 

To get the air mass flow (AMF) from the volumetric flow, an estimate of 
air density is needed. Using the assumptions that the air was relatively 
dry (the relative humidity was under 50%) and behaves as an ideal gas, the air 
density was calculated as: 


Density = 


Pressure 

air gas x absolute 
constant temperature 


The plenum pressure was small (.02 psi) so the absolute pressure was 
2120 pounds per square foot — close to atmospheric. In English units 


. 2120 3 

Density = pounds mass/ft . 

53.3 [ 460 + T plen „ <®F)J 


Since air mass flow is density times volumetric flew, 
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52.224 x 10 3 

AMF pounds mass per minute. 

^ + T plenum ^ 

Over the range of temperatures encountered in the input plenum, which was 
59°F to 71°F, AMF ranged from 100.62 down to 98 .35 pounds mass per minute. 
AMF is shown in the data sheets. 

The thermal power output of the heat pump, either in the heating or cooling 
mode, is: 

P QUt (thermal) = AMF x specific heat x temperature difference. 

The appropriate specific heat constant is that for dry air at constant 
pressure since the pressure change in the furnace is very small. The value 
is Cp = 0.24 BTiyifcm-°F so the output formula becomes: 

P out (thermal) = AMF x C p x AT(°F) = 14.4 x AMF x AT(°F) BTq/hour, 

which is in the units that heat pump manufacturers conventionally use. When 
this BTt/hr figure is divided by 3413 BTq/hour, the power output is in units 
of kilowatts. When the power output in kW is divided by the electrical 
power input in kW, the resulting number is the coefficient of performance 
which is about 3 for the test conditions. The power output calculations at 
nominal conditions yielded power outputs close to the manufacturers specifi- 
cations. 


Heat pump io coding mode: 

In this mode, where only a few measurement runs were made, the furnace 
fan was connected to the variable source so the air flow thru the furnace 
was affected by the source. Over the frequency range of 50 to 70 Hz, the 
fan plenum, pressure changed only about J. inches of water or less than one 
part in 10 3 in absolute pressure so the air mass flow will change 
significantly only due to volumetric flow changes. If one assumes this flow 
change will be proportional to frequency only, then 


Fan frequency 

Air flow = Air flow (60 Hz) x -- 


Two checks were made at 70 Hz to validate this assumption. The measured 
peak velocity out of the exhaust vent increased to 2500 ft/min from the 2200 
peak value at 60 Hz. If the average flow increases as much, the correction 
2500 is 1.14 while the frequency correction just described is 70/60 or 1J.7. 

2200 

The second check used the increased pressure drop of 0.40 inches of water 
across the "A" coil at 70 Hz where it was 0.31 at 60 Hz. This suggests a 
correction factor of 


28 



1/2 


AMF (70 Hz) = AMF (60 Hz) 


x 


AP (70 Hz) 
AP (60 Hz) 


= AMF (60 Hz) x 



= AMF (60) x 1.136 


The square root is used since AMF is related to velocity (at constant 
density — already justified) and velocity change is proportional to the 
square root of pressure change from the energy equation for compressible 
flow. The checks are OK so the output power of the heat pump in cooling 
mode with furnace fan excited by the variable frequency source is 

P out (thermal) = 14.4 x AMF C x AT(°F) BTCJ/hour 

where AMF C is AMF corrected by frequency ratio. The AM^and power output data 
obtained by the calculations described above are included in the data sheets 
(Appendix A) . 
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V. EXPERIMENTAL RESULTS 


Since data was available for about 170 test runs it would be easy to 
overload the reader with multiple presentations of data. Furthermore, the 
conventional method of presenting motor data [eg: Ref. 2] , is to plot with 
motor shaft power output as independent variable — and the motor shafts were 
inaccessible for measurement.' Thus, the investigators were forced to pre- 
sent the data in an unconventional manner and, after considerable delibera- 
tion, decided upon the "operating envelope" concept common in aircraft 
performance analysis. Here the range of variables for stable operation is 
presented with the suggestion that operation outside of this envelope is 
either unstable or unsafe. The usual variables for aircraft in level flight 
are speed and altitude while in this report, the investigators chose input 
voltage and frequency. In some plots, the power input and power output 
corresponding to each voltage- frequency data pair are included. A glance at 
the "envelope" is instructive as it suggests the range of input combinations 
at which the particular pump will operate. It is also possible to indicate 
the reason for unsuccessful operation at the edges of the envelope. The 
reader will notice that motor stall due to inadequate torque, in turn due to 
inadequate excitation voltage, is noted at the bottom of the envelope and 
cavitation and overload is often present at high frequencies. High voltages 
would be expected to also increase iron losses with resultant overheating. 
However, most runs were too short for overheating to appear and besides, 
winding temperatures were not measured. 

It should be noted that the high voltage side of the envelopes are not 
completely explored since the excitation system of the source, even with the 
autotransformer is inadequate to provide voltages much beyond 50% of nomi- 
nal. The other presentation technique used was a conventional plot of load 
power ratio (where nominal or rated is unity) as either voltage or frequency 
is varied while the other is held nearly constant. All of the measured data 
points come from the data sheets in Appendix A. 


l) Experimental results for Myers Z 1Q-1§ submersible pump 

In Fig. 10, a partial envelope of the submersible pump indicates that 
this pump will operate between 42 and 100 Hz and requires less voltage as 
the torque requirements decrease with pump speed (frequency). The power 
input and output numbers show that the overall efficiency decreases above 75 
Hz when apparently cavitation begins. There are too many data points 
closely grouped to include the power numbers. Therefore, just the voltage 
and frequency pairs are plotted in Fig. 11. A stable condition means that 
the measured pump variables did not change significantly during that test. 

In Fig. 12, the submersible pump power ratio from nominal frequency is 
plotted for voltage within 3% of nominal 230 volts. A linear representation 
of this data indicates a change of pump power input of 2.2% for a 1% change 
in frequency. The increase in power with cccfrequency indicates that as 
windmill rotor speed increases during a gust, the load will also increase to 
moderate rotor acceleration. 


30 




Measured. Partial Operating Envelope of Submersible Pump 
Shotting Input and Output Power 
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Measured Stable Voltage-Frequency Pairs for Submersible Pump 
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Power Ratio 



Figure 12 


Variation in Submersible Pump Power 
as Frequency is Changed 


2) ExpetimeoUl results fee E*E t Myets QP1Q. ceotcifugel citculatiog eu®p 


The measured operating envelope for the circulating pump with long pipe 
load is shown in Fig. 13. The frequency range is less than the submersible 
pump— from 48 to 92 Hz but the stable voltage range is greater — from 94 to 275 
volts. The same pump with a lesser load (pipe length about 1/2 of above, gate 
valve wide open and no inlet valve) has a higher output but operates over the 
same range indicating the range is dominated by motor characteristics. Not 
enough data points are available for a reasonable estimate of the operating 
envelope under this loading. The data available is shown in Fig. 14. The two 
figures indicate that the centrifugal pump can tolerate more voltage variation 
than the submersible pump but that stable operation at low frequencies is 
harder to maintain. Furthermore, this type of pump requires a foot valve at 
its suction pipe inlet to avoid loss of prime. The investigators suspect that 
maintaining prime (avoiding loss of suction) will be required for intermittent 
operation as provided by a windmill source. Fig. 15 illustrates the behavior 
of this pump as voltage is changed while the frequency is maintained within 
1/2 Hz of 60. The pump power input is not sensitive to voltage changes so 
that 1% increase in voltage causes only a 0.4% increase in pump input power. 


3) Experimental results far Addison Q0&1& beet pump 

Unlike the water pumps, the heat pump was field tested under reasonable 
climatic conditions with a test setup capable of greater excursions in voltage 
and frequency and with greater redundancy in output measurements. Fig. 16 
shows the measured stable operating points for the heat pump in heating and 
cooling modes. The addition of the autotrans former, not available for the 
water pump tests, allowed voltages to reach almost 400 volts at the higher 
frequencies and the pump operated without noticeable complaint down to 50 Hz 
before stalling around 47 Hz. The lowest voltage of operation is 175 volts 
which confirmed a test run by the Principal Investigator in June 1971 when 180 
volts was the lower limit of stable operation for a Sears 14,000 BTt/hour air 
conditioner rated at 230 volts [Ref. 8, p. 29]. In that 1971 test, the 
compressor- motor power factor was close to unity at that voltage which is also 
true for the heat pump at low voltage. The 1971 tests were run only at 60 Hz. 
The Addison chief engineer warned that the compressor outlet temperature must 
remain under 260° F, to avoid lubrication problems, and the data indicates a 
maximum of 213 °F which occured during a 90 Hz run. Fig. 17 shows a subset of 
the data with input and output numbers included. Note that the coefficient of 
performance approaches four near 50 Hz but drops toward two over 85 Hz. and at 
voltages over 320 where motor excitation and pumping losses increase. When 
the voltage decreases, the motor current increases so winding losses increase 
and high motor temperatures seem likely. This was verified by the opening of 
the thermal cutout in the motor winding which occurred after ten minutes of 
operation at low voltages. The motor was drawing (both the compressor motor 
and outside fan motor together) 13 J. amps (rated is 11 amp6) when the cutout 

opened. Nevertheless, the heat pump will operate on a short-term basis over a 

range of voltage and frequency that surprised the investigators. 

% 

The variation in heat pump input and output power as frequency is varied 

is shown in Fig. 18. Note that the power output line is much flatter than the 

input line and the crossover is at 92% of rated frequency which support the 
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Figure 14 

Measured Partial Envelope of Circulating Pump uith Shorter Pipe Load 
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Figure IS 

Variation of Circulating Pwnp Input 
Power as Voltage is Changed 
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The variation in heat pump input and output power as frequency is 
varied is shown in Fig. 18. Note that the power output line is much flatter 
than the input line and the crossover is at 92% of rated frequency which 
support the conclusions made earlier in this section on the variation of 
coefficient of performance. The variation of heat pump power with voltage 
is shown in Fig. 19 but a reasonable linear fit of the data is unlikely. 

The power input data for 60 and 80 Hz have the same shape as in the 1971 
tests referenced earlier. The only other test data to report is summarized 
in Fig. 20 where the current waveforms of the heat pump are sketched for 
several extreme and the nominal condition. Note the magnetic circuit 
saturation characteristic at high voltage and frequency and the triangular 
waveform due to the odd harmonics caused by excitation current at lower 
voltages and frequencies. These sketches were made by Dr. J. Rumba ugh, 
Program Manager for the DOE Wind Energy Technology Division during his site 
visit the last day of testing. May 21, 1985. 
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VI. SIMULATION METHODOLOGY 


The simulation approach chosen represents the single- phase caged motors 
with the classic revolving field equivalent circuit, including the start- run 
capacitor, as shown in Fig. 23. The pump load is modeled using the 
mechanical power and slip (related to motor shaft speed) variables present 
in the equivalent circuit. The resulting electrical and mechanical 
equations are then programmed on a small computer in PASCAL. The output 
plots from the computer runs are shown and explained in the next section. 


The investigators did not have time to obtain the equivalent circuit 
data for the heat pump motors as well as the model of the refrigerant loop. 
Furthermore, not enough measurements were available to characterize that 
loop. Therefore, no simulation was made of the heat pump behavior. 
Fortunately, the input-output measurements on the heat pump were extensive 
so simulation was of lesser importance than with the water pumps. 

1) Pump Hydcsulic Loads 


The experimental setups for the circulating and submersible pump tests 
are shown in Figures 3 and 1 respectively. Assumptions of wide open gate 
valves and neglect of the loss across the flow meter and the gate valve were 
made. The standard formulation for friction head loss in a pipe of length 
(ft), inner diameter D (ft) with flew Q (f ty sec) is 



Q 2 

2g(D 2 ty'4) 2 


(ft) 


where g is the gravitational constant and f is the Moody friction factor. 
The following closed form approximation [Ref. 10, p. 292] for f is 

f = a + bR e " c 


where 


= the Reynolds number - 

® Tll/D 

v = the kinematic viscosity 
a = 0.094 k 0,225 + 0.53 k 
b = 88 k 0 * 44 
c = 1.62 k 0 * 134 

and 

k = relative surface roughness, assumed equal to 10“ 5 for plastic 
pipe. 
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Hie circulating pump setup has a static head of 6 feet and contains 60 ft 
of 1 1/4" ID plastic pipe and 345 ft. of 1" ID plastic pipe. The submersible 
pump setup has a static head of 28 ft. and contains 252 ft. of 1" ID plastic 
pipe. The sum of static and friction head theoretical loads on the two 
pumping setups are shown in Fig. 21. 

2) Pump Models 

The standard "pump curves" supplied by the manufacturer for the 
circulating and submersible pumps are shown in Fig. 22. In actuality, these 
are not truly "standard" curves as commonly used by pump manufacturers as 
they are oot at a constant speed. Rather, they are as driven by a given 
induction motor which only has approximately constant speed operation. For 
simplicity assume that these curves ace at a constant speed, the nominal 
operating speed at rated output of the motors, 3500 rpm for the circulating 
pumps and 3440 rpm for the submersible pump. These pumping curves were 
approximated by the closed form expression 



for the circulating pump and 



for the submersible pump. In these expressions, Q is in units of gpm. Hie 
pump efficiency curves were approximated by the closed form expression: 

eff c = .658 sin (j £--) 

for the circulating pump and 

eff s = .518 sin (j 

for the submersible pump. Again Q is in gpm. 

For speeds off nominal, the hydraulic affinity relationships were used 
to find H vs Q for points of equal efficiency to those of the equations 
above for H c and H s . 

That is 

Q= _^m 

^ ^mo 
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Figure 2.1 

Theoretical Hydraulic Pumping Loads 
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Figure 22 

Pump Curves for the Water Pumps 
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and 


2 



where f m is pump speed and the subscript o refers to head and flow at the 
standard speed f^. 

3) Eump {totQc Models 

Manufacturer's data for the permanent split-phase capacitor (PSC) 
squirrel cage induction motors supplied with the pumps used in the tests 
are shown in Table 2. These machines are modeled with the double revolving 
field theory as detailed in [Ref. 11]. The revolving field equivalent 
circuit derived for a PSC motor is shown in Figure 23. The input voltage 
sources for the forward and backward revolving fields are given by 

V f = V s (l - j/a)/2 


and 


v b = V s (l + j/a)/2 

where a is the effective turns ratio of the stator capacitor winding to the 
stator main winding and V s is the applied RMS voltage to the electrically 
paralleled stator windings. The total stator input current is the phasor sum 
of the two stator winding forward and backward currents transformed back to 
their respective stationary components, and is given by 

z s * I fs (1 + i/ a) + I bs {1 “ i /a) * 

Thus, if V g is the phasor reference, the input power to the motor is 

p in “ v s x ^al (I s ) , 


Real (I s ) 

the input pcwer factor is PF = 


and the mechanical output power is the difference between the forward and 
backward powers imparted to the rotor, that is 

2(11*1**!=!- U br l 2 ) 

where the factor of two is due to the definition of the revolving field 
transformation. 

Thble 5 is the list of equivalent circuit element values for the two 
water pump motors as determined by "trial and error" fitting of element 
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The upper loop is the forward revolving field circuit and 
the lower loop is the backward circuit 


R 

r 


1-8 


8 


R I=£ 
r 2-s 


The parameters are: 

Z , - %( Z /a 2 - R - joj/ ). and Z = R + jai L + -t~ 

a as s s as as ^ as jii £ 


as 


a = ratio of effective main stator winding turns to 
capacitor winding turns 

R = stator main winding resistance 
s 

h as = stator capacitor winding resistance 

L= stator main winding leakage inductance 

T as = stator capacitor winding leakage inductance 

C = start-run capacitance 
as r 

L^= magnetizing inductance (RTS) 

R = equivalent resistance representing magnetizing losses 
w (pj’g) 

L^= rotor leakage inductance (RTS) 

R = rotor resistance (RTS) 
r 

s = rotor slip = (f - f m Vf 

f = electrical supply frequency in per unit (60 Hz base) 

/ = mechanical speed in per unit (3600 rpm base) 

(RTS) means referred to stator 


Figure 23 

Equivalent Circuit for PSC Caged Motors 
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Element 

Chit 

Circulating 
Pump Motor 

Submersible 
Pump Motor 

a 

numeric 

1.25 

1.25 

*s 

-tl 

2.43 

1.814 

Kas 

-TV 

4.0 

3.0 

U>o L s 

JO- 

2.5 

3.0 

^o L as 


3.0 

3.5 


.O- 

106 

106 



45.0 

37.0 



500 

225 

Wo L r 


2.0 

2.0 

«r 

JTL~ 

3.8 

3.5 

CO o 

rad/sec 

377 

377 


able 5 

Circulating Pump Motor and Submersible Pumps 
Motor Equivalent Circuit Element Values 
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values such that simulated performance curves matched those provided by the 
motor manufacturer or matched test points so provided. The equivalent 
circuit values in Table 5 yielded the pump motor characteristics shown in 
Figs. 24 and 25 where motor efficiency, speed and power factor are plotted 
in % against motor shaft power P m and where motor current is in amperes. 


The circuit model for motors is assumed linear so the magnetizing 
inductance L_ is assumed constant. As a consequence, the motor output P m 
has a square-law dependence on the magnitude V s of the stator voltage V g 
given by, 

^m^s' = v s *ml^ f ^ ' where fnj) 


is the mechanical output power of the motor as a function of the electrical 
frequency f and the shaft speed f m , where the latter is taken at nominal input 
voltage of 230. V s is in per unit (percentage of 230 volts divided by 100). 

For a hydraulic load of fixed configuration, the mechanical shaft speed f m is 
the only independent variable in the pump and load equations. Thus it is 
possible to display contours of constant pump output or flow i.e.: constant 

pump speed due to the fixed hydraulic load configuration, on the same stator 
voltage- frequency plane chosen to display the operating envelope. Here 
specifying f m determines pump output and pump efficiency and, indirectly, the 
stator voltage V g via the equation 


’s 


Pump output (f^ 
efficiency (f m )' (f, f^j 



p.u. 


These contours are exhibited in the figures in the next section. 


The pump output Ppumo on ^ a function °f the pump speed f m but it 
must be determined in xorrjunction with the hydraulic load. Once the speed 
fjp has been specified, the pump head versus flow relationship is solved 
simultaneously with the load head versus flow relationship to find their 
intersection or operating point. This process for the sumbersible pump 
setup is shown in Fig. 26 where the intersection of the load curve with pump 
curve determines the total head and flow for a particular pump speed, which in this 
case, is the rated pump speed f^. 

The sequence of calculations to generate the voltage- frequency contours 
indexed by flew is then as follows: 


1) Specify the mechanical speed f m . 

2) Generate a new H vs Q "pump curve" for the speed f m by use of the 
hydraulic affinity relationships. 

3) Simultaneously solve the H vs Q "pump curve" and the H vs Q load curve 
for an operating point H and Q. From this point calculate the pump 
output power. 
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4) Determine the operating point pump efficiency from the rated speed pump 
efficiency curve and the hydraulic affinity relationship between flow and 
speed. 

5) Calculate, as a function of electrical frequency f, the required motor 
stator voltage from the square root relationship between V g and the pump 
output, pump efficiency and the motor power for rated (1 pu) stator 
voltage. 

6) Plot the result in the V g , f plane 
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VII. SIMULATION RESULTS 


Calculated results for constant pump output contours in the pump voltage 
frequency plane are given in Figure 27. The smooth contours are labeled both 
in pump speed f m and pump flow Q (gpm). As can be seen, for any speed f m , 
steady-state operation is obtainable only for electrical frequencies f > f m 
(both f and f m are in per unit), that is, for positive rotor slip. As f 
approaches f_, the stator voltage V g must rise as the equivalent rotor 
resistances also rise with decreasing slip. Conversely, as f grows above f m 
the stator voltage falls, since the equivalent rotor resistances also fall 
with rising slipw The contours are stopped when the slip reaches 20% since 
for these particular motors, their torque output is maximum at that slip. 

Also shown is Figure 27, are the four experimental points for which the 
gate valve was wide open. Correlation between theory and experiment for these 
four points is quite good with only the high frequency point being off by more 
than one gpm. 

A similar set of theoretical constant pump output/speed contours for the 
circulating pump setup is shown in Figure 28. Here again the contours are 
halted at the value of electrical frequency f which gives the "pullout" slip 
of 20%. However, in the circulating setup there were a greater number of 
suitable, i.e., open gate valve, experimental points to compare to the 
theoretical prediction. These experimental points are also shown in Figure 
28. Again, the agreement between theory and experiment is quite remarkable 
and again the greatest discrepancy is at the high frequency point. 

Also shown in Figure 28 are two straight line contours of desired 
operation. One, labeled "nominal rated operation," is the desired operating 
locus if one had no information on the details of the operating system and the 
other, labeled "Max. eff. operation," is the optimum operating locus for the 
given fixed hydraulic load, specific pump and specific motor. This 
difference in operating strategies is more obvious in Figure 29 where the 
circulating pump motor efficiency is shown as a function of mechanical speed 
f m for different electrical frequencies f. Recall that the motor model is 
assumed linear and as such the motor efficiency is not a function of motor 
terminal voltage. Operation at the maximum efficiency point for a given 
electrical frequency is then only a question of operating at a specific motor 
speed or slip. These optimum slip values are those above the maximum 
efficiency contour in Figure 28. This contour of optimum slips then also 
determines the optimum voltage divided by frequency, or vVfr contour for the 
specific pump load. If one were to operate at the nominal V^/f ratio of unity 
(both V s and f in per unit) as shown in Figure 28, the operating slip at the 
different values of speed can be extracted and then replotted on the efficiency 
curves of Figure 29. This contour clearly shows that operation at a nominal 
Vg/f ratio of unity is inferior, particularly so at frequencies below rated. 

If one assumes an optimum wind alternator field controller, one that keeps 
the (V , f) steady-state operating contour along that shown in Figure 28, the 
max. erf. contour — the minimum needed pump motor input power as a function of 
input frequency — can be uniquely determined. This calculated input power, for 
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Figure 27 

Simulated Voltage-Frequency Contours 
for Submersible Pimp and Load 
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Figure 28 

Simulated Voltage-Frequency Contours 
for Circulating Pump with Load 
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Figure 29 


Mechanical Speed, f m (pu) 


Simulated Efficiency vs Shaft Frequency 
for the Motor in the Circulating 
Pump 
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6 




the circulating pump motor , is shown in Figure 30. Note that this minimum 
input power contour is a close fit to a cubic relationship in operating 
frequency. Since the power available from the wind is also a cubic function 
of the wind speed for a given swept area, this load is ideal to maintain a 
constant wind machine power coefficient C_ (dependent only on the ratio of 
wind speed to wind machine rotational speed, i.e., f for an alternator). 

Thus, a properly controlled alternator, induction motor, and centrifugal fluid 
pump working a fixed hydraulic load is an ideal wind machine load, one that 
enables maximum steady-state power extraction from the wind over a fixed wind 
speed range [Refs. 5,12]. 
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VIII. REQUIREMENTS FOR WIND DRIVEN GENERATORS AND CONTROLLERS 

FOR ISOLAIED OPERATION 


To supply varying electrical loads with constant frequency and voltage 
requires generation capacity larger than the maximum load and using precise 
governors and a voltage control feedback system; usually including exciter 
feedback control, automatic tap-changing transformers and switched capacitor 
banks. A well-designed and operated system of several generators and many 
varying loads usually maintains frequency within 1% of 60 Hz and the service 
voltage within 5% of nominal unless there is a component failure or unless a 
relatively large varying load is connected. When interconnected with large 
utility networks, the frequency variation of a small utility as just des- 
cribed will significantly decrease but the voltage changes will still be 
approximately 5% or less. When a small controllable source, such as an 
isolated single generator driven by an internal-combustion engine is used, 
it is difficult to keep the voltage and frequency variations small unless 
the electrical load variables are constant. Nevertheless, a small auxilary 
alternator with a good governor and voltage regulator will keep voltage and 
frequency within the normal limits specified by most electric motor manufac- 
turers. However, when the source of alternator power is itself highly 
variable, such as from a wind-electric system, then the problem of control- 
ling frequency and voltage within tight limits becomes impossible under 
isolated conditions. Even when an electric windmill is interconnected to a 
strong utility, so that the frequency is tightly regulated, the fluctuations 
in wind velocity cause corresponding, but much smaller, fluctuations in 
local voltage while the power injected by the windmill has peak fluctuations 
often equal to its 15 minute average power [Ref. 9, Fig. 7], Thus it seems 
reasonable to state that supplying a varying load with an isolated windmill 
without battery storage will result in extensive frequency and voltage 
variations even when sophisticated controllers and/or power conditioners are 
used. Since many self-excited wind-driven alternators will operate over a 
wind speed range of at least four to one under stand-alone or isolated 
conditions, it seems unlikely that shaft speed (and thus frequency) and 
output voltage can be held within 5% of a nominal speed without unacceptable 
performance compromises. This conclusion is of course the premise of this 
study and the reason for the tests to estimate the performance envelope of 
critical agricultural loads. 

It is apparent from Sections V and VII that the loads considered in 
this study will only operate over a frequency range of approximately 50 to 
90 Hz and a voltage range less than 100 to 350 volts depending upon the 
mechanical load characteristic. Thus, since the wind speed variations 
normally encountered are likely to generate voltages and frequencies that 
exceed the observed limits, controllers will be necessary on windmills used 
to feed such loads. It is important to note here that the common induction 
generator used in many utility-interconnected windmills requires essentially 
no operating controls except for a contactor. The purpose of this section 
is not to specify a controller for an isolated wind-driven alternator but to 
note some of properties of such a controller and its associated windmill 
suggested by the tests and analysis performed in this study. 
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Properties needed fee voltage control r 


The magnitude of the voltage across the load, which is assumed to be 
not far from the windmill so transmission effects can be ignored, is 
proportional to alternator shaft speed, magnetic flux due to field current 
and the turns ratio of any transformer used between alternator and load. 

Capacitive shunt compensation across the load can also increase voltage 
significantly if alternator stator inductance is large. The envelope plots 
in Figs. 10, 11, 13 and 14 indicate that the water pumps will operate with 
much less than nominal voltage at lower frequencies, although their output 
is higher if voltage is kept at a voltage to frequency ratio (V/f) equal to 
that at rated conditions (V/f = 3.8). For all devices tested above rated frequency, 
the control of voltage so thatV/f is held at nominal, does not increase the output 
significantly and leads to excessive losses and thus heating in the motor 
windings. This suggests a first requirement: 

* The alternator field circuit (including magnetics), voltage 
regulator and any auto-switched transformers should be designed 
so the voltage is reduced as rotor speed reduces so V/f is 
constant below rated frequency and V is near rated above 
rated frequency. 

Furthermore, a high and low voltage cutout function should be included if 
V/f control is not implemented to avoid operation at voltages below rated V/f 
where stalling occurs and much over rated voltage at high frequencies where 
motor overheating occurs. This suggests: 

•Implement a high limit 20 to 30% over rated voltage and a low limit at 

the voltage below rated where stall due to lack of motor torque occurs. 

At low frequencies, motor torque decreases due to decreased capacitor- 
winding current which could be compensated by increasing motor capacitance 
at lew frequencies. 


Properties needed foe frequency control! 

Depending upon the type of load and any capacitor switching, the load 
must be disconnected between 40 and 50 Hz to avoid stall and subsequent 
overheating. To avoid cavitation and pump overload, operation above 90 Hz is 
undesirable. Since normal wind speed variations at most sites will result in 
larger excursions in frequency, either blade pitch control or resistive load- 
ing and unloading is necessary. A resistance (heating) load could be switched 
in at high frequencies and switched out at low frequencies to help control 
the shaft speed variations. 

* Implement either pitch control or controllable resistance loading to 
reduce the frequency range to about 50 to 90 Hz. Switch out pump loads 
at both high and low frequency limits and add resistive load at high 
frequencies to avoid rotor overspeed when the pump loads are switched out. 
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Fortunately, the water pump load power characteristics used increase with 
rotor speed in a cubic manner similar to that of rotor shaft power output so 
centrifugal water pump loads should not require much control within the fre- 
quency limits discussed. 

Etopecties needed foe pump starting; 

Motor starting currents (locked rotor amperes) are 6 to 8 times rated 
currents due to motor electrical behavior and to the high inertia of the water 
mass or refrigerant pressure in which the pumps are immersed. There is no 
practical inexpensive way to unload the pumps at start without extensive 
modifications of the pumps as well as loss of prime and lubrication, so high 
starting currents must be provided by the alternator. Implementation of this 
function can be complex. One required electrical property of the alternator 
is: 


* The alternator synchronous impedance (Thevenin equivalent impedance) 
must be low enough so the alternator short circuit current exceeds the 
locked rotor current of the loaded pump motor. 

This requirement may prove challenging when using a large diameter many- 
pole alternator equal to the pump power rating which will generate 50 to 90 Hz 
without a speed-up gear train. Since the loads approach operating speed 
within a second or two, the possibility of starting with a wind gust appears 
attractive. This suggests a high- inertia wind rotor with a speed transducer 
and a control algorithm that waits until enough kinetic energy is accumulated 
in the rotor to start the load and still have enough speed to stay within the 
operating envelope of the load. 

* Implement a load starting controller that monitors rotor kinetic energy 
until a start with subsequent operation is likely. An algorithm adapted 
to the local wind regime is likely needed. 

Miscellaneous controller properties suggested to the tests i 

* Add thermistors to motor or pump casings to disconnect loads at housing 
temperatures that are likely to open thermal cutouts, reduce motor life 
significantly, or adversely affect refrigerant- lubricant properties. 

* Add an appropriate time interval to the cutout algorithm so the loads 
can reestablish equilibrium of the pumped fluid or gas. 

Additional studios needed: 

Although several pumping loads were successfully operated far outside 
normal voltage and frequency limits, there are several investigations neces- 
sary before the approach in this report can be recommended. They are sug- 
gested by the limitations adopted for this study: 
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Transient operation with windmill source. 

The operating envelopes were determined from several-minute steady-state 
runs with the source at constant frequency and voltage. Using wind 
generation in a real wind regime will require operation of critical loads 
and source under transient conditions including starting and possible 
wind-turbine stall. 

Use of loads driven by motors with switched starting windings. 

Long-term studies estimating effects of continued large fluctuations in 
frequency and voltage on motor and load service life and reliability. 
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Appendix A 

Transcription of Field Data 
Sheets Including Calculated Values 
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CThis Turbo Pascal program will perform the analysis for the centrifugal 
pump psc motor, data given m Table 5, and plot on the screen the results 
shown in figure 25. Turbo Pascal is a trademark of Borland International , Inc . } 


CThe following pascal functions perform network analysis. They are roughly 
based an the ane-/two-port linear circuit analysis program MARTHA} 

type 

vector - array C 1 .. 2000} of real; Cratating memory for complex quantities} 


var 

f 

rl,r2,r3,r4 

r 

col , cc2 
kkl , kk2 


real ; 
integer ; 
vector ; 
integer ; 
integer; 


const 

dumi : integer * 1 
dumiS : integer - 3 
dumi3 : integer - 5 


function cl : integer; Cpointer to the real part of a complex number} 
begin 

ccl:-ccl+2: if Cccl>-1993 then ccl :- 7; 

Cl CCl'; 

end; 

function c2 : integer; Cpointer to the real part of A of the ABCD parameters 

of a two-port network} 

begin 

cc2:“Oc2+B; if Ccc2>-9923 then cc2 : ” 201; 
c2 ccB ; 

end; 

function kl : integer; Csimilar to cl but designated keeper space} 
begin 

kkl: -kkl +2; 
if Ckkl>-11993 then 
begin 

writelnC ’keeper counter kl has rolled over’}; 
kkl: -1001; 
end ; 

kl :- kkl; 
end; 


function k2 : integer; Csimilar to c2 but designated keeper space} 
begin 

kk2: -kk2+8; 
if Ckk2>-19923 then 
begin 

writelnC 'keeper counter k2 has rolled over’}; 
kk2 : -1201 ; 
end ; 

k2 :- kk2; 
end; 

function caCi,J,k : integer! : integer; Ccomplex addition} 
begin 
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rCIO : “rC i!*rC J ! ; rCk*13 : -rCi+ll+rC J*13 ; 
co : " k; 
end; 

function csCi,J,k : integer} : integer; Ccomplex subtraction! 
begin 

rCk3:-rCi3-rCj3; rCk*13 : -rCi+13-rC J+13 ; 

CS : “ k ; 

end; 

function cmCl,J,k : integer! : integer; Ccomplex rrultipl icotion) 
begin 

rCk3:-rCi3TCj3-rCi+13“rCj*13; 
rCk+13 : -rCi3»rCj*n*rCi+ll»rCj3j 
cm : “ k ; 
end; 


function cdCi,J,k : integer! : integer; Ccomplex division! 
var dum ; real; 
begin 

dum;-rCj!»rC J3-TC j*13“rC J+13 ; 

rCk! : “CrC i!*rC jl+rC i*l!*rC j+1! 1/dum ; 

rCk*13: -CrCi<-13“rCj3-rCi3»rCj-'-13 3/dum; 

Cd : “ k ; 

end; 


function ciCi.J : integer! : integer; Ccomplex inverse! 
var dum : real; 
begin 

dum: -rC i!*rC i!+rC i*13*rC i*13 ; 

rC j! : - rC i! /dum; rC J+13 : —rC i*13 /dum; - 

Ci J; 

end; 

function ccCi.J : integer! : integer; Ccomplex conjugate! 
begin 

rCj3:-rCi3; rC j+13 ; — rCi+13 ; 

CC J; 

end; 

function cmplxCx.y : real; i : integer! : integer; Cforms a complex number! 
begin 

rCl!:-X; rCi*l!:“y; 

cmplx i; 
end; 


function zrCr : real! : integer; Cforms a resistive one-port! 
begin 

zr cmplxCr,0.0,cl!; 

end; 

function zlCl : real! : integer; Cforms an inductive one-port! 
begin 

zl cmplxCO .0, f *1 , cl! ; 

end; 

function zcCc : real! : integer; Cforms a capacitive one-port! 
begin 

zc cmplxCO . 0, -1 .0/Cf*c! , cl! ; 

end; 
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function sCi,J : integer] : integer; { series connection of two one-ports ] 
begin 

S: -coCi , J , Cl] ; 
end ; 

function pCi.J : integer] : integer; {parallel connection af two one-ports] 
begin 

p: -cdCcmCi , J,dumi!,caCi,J, dumi23 , cl] ; 
end ; 

function wsCi : integer] : integer; {creates a tuio-port from a series 

connected one-port! 

var d : integer; 
begin 
d:-C 2; 

rCd! : ”1 .0; rCd*63:-1.0; 

rCd*l! : -0 .0; rCd*4] : “0 . 0; rCd+5! : -0 . 0; rCd+7] : -0.0; 
rCd*2! : “rC i! ; rCd*33 : -rC i+13 ; 

US d; 

end; 

function wpCi : integer] : integer; {creates a two-port from a shunt 

connected one-port] 

var d : integer; 
begin 
d: “C2; 

rCd] : “1 .0; rCd+S] : -1 . 0; 

rCd+1! : ”0 .0; rCd*23:-0.0; rCd+3!:-0.0; rCd+7] : -0 .0 ; 
up: “CiC i , d + 4] ; up d; 
end; 

function wcli.J : integer] : integer; {cascades tuo two-parts! 

" var d : integer ; 
begin 
d : -c2 ; 

d: -calcmCi , dumi] , cm{i+2, j+4, dumi23 , d] ; 
uc:“caCcmti, J*2,dutni! , cmC i+2 , j*6 ,dumi2] , d+23; 
uc: -calcmCi+4, J , dumi] , cmli+6, J+4, dumi2! , d+4] ; 
uc : -caCcmC i*4 , J+2 , dumi ] , cm{ i+6 , J*6, dumi2] , d*S] ; 
uc d; 
end; 


function utsli : integer] : integer; {terminates a two-port with a short! 
begin; 

uts cdC i*2 , i+6, cl! ; 

end; 

function wtoCi : integer] : integer; {terminates a two-port with an open] 
begin 

wto cdCi, i+4,cl3 ; 

end; 

function wtCi.j : integer! : integer; {terminates a two-port with a one-port! 
var d,dd : integer; 
begin 
dd : -cl ; 

d : ”caCcm{ i , j , dumi ] , i*2 , dumi2] ; 
ut cdCd , caCcmC i<-4 , J , dumi ] , i+S , dumi3! , dd] ; 

end; 
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function kesplCi : integerl : integer; Cputs a complex number in keeper space! 
var d : integer; 
begin 
d;-kl; 

rCd3;“rCil; rCd«-l! ; -rC i*l! ; 
keepl d; 
end; 


function keep2Ci : integer! integer; Cputs o two-part's PBCD parameters 

in keeper space! 

var d,dd : integer; 
begin 
d:-k2; 

for dd;“0 to 7 do rCd*dd! : - rCi*ddl; 
keep2 d; 
end; 

procedure c_initial ize ; Cinitialize normal pointers! 
begin 

CCl:“7; CC2 : ”201 ; 
end; 

procedure k_initialize; Cinitialize keeper space pointers! 
begin 

kkl:-1001; kk2 : “1201 ; 
end; 


CThe following are Turbo Pascal screen graphics/plotting functions 
for an IBrt-PC with a graphics adapter.! 

var - - 

i , j : integer; 
sx,sy :real; 

procedure axies; 
begin 
hires ; 

hirescolarCyellowl ; 
for i : “5 to 505 do plot C i , ISO , 1 1 ; 
for i;“l to 190 do platC5,i,13; 
end ; 


procedure ticksl; 
begin 

for i:“l to 10 do 
begin 

J:-5*B0«i; 

platC J, 189, 1! ; plotCj, IBB, 11; 
J : “190-i*19; 

platCB, J, 11; plot C7 , J , 11; 
end; 
end; 

procedure gridl; 
begin 

for i:“l to 10 do 
begin 
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* 


For J:-l to 600 do plotC5*J , 190-19»i , 1 3 ; 
For J:-l to 190 do plot C5+i»60 , J , 1 3 ; 
end j 

end; 


CThe Following is the main poscol program to analyze the psc motor . 3 
var 

rs, ls.ras, los.cas, lm, lr,rr,rm 
zrb , zrF , zmb , zmF , zd , zsF , zsb 
k , m, n, x, y 

slip , Fm, zbase, vbase , pbase , pm 
yfF,yfb,ybb,ybF, vf , vb 
a , aa , vs , pF , pb , mirF , mirb 
zm, zdF , zdb, zs , irF , irb , zo 
isF , isb, denom, seFin, sebin 
pein , eFF , mi_line , power_Factor 
im, ia, i_line 

begin 

Cmotar constants! 

F:-1.0; rm:-500.0; 
rs:“2.H3; rr:-3.B; 

Is ; “2 . 5 ; lr:-2.0; 
lfTU-HS.O; cas:-l. 0/106.0; 
ros^'i.O; las:“3.0; 
a : “1 . 2S ; aa:-2.0*a*a; 

VS: -230.0; 

axles,- gridl; 
pbase : -2 . 5»7 l i6 . 0 ; 

SX : —600 ; sy. “190; 


real ; 

integer ; 
integer ; 
real ; 
integer ; 
real ; 
integer; 
integer ; 
real ; 
inteoer : 


For i:- 1 to 150 do * " 

begin 

Fm : -C 1000-i 3/1000 ; 
c_initialize; k_initialize; 
slip: -CF-Fm3/F ; 

Csetup the network shown in Figure 23.3 

zrb : -sCzlC lr 3 , zrCrr/C2 .0-slip 3 3 3 ; 

zrF:-sCzlClr3,zrCrr/slip33; 

zm: -pCzrCrm3 , zl C lm3 3 ; 

zmb: -pCzm,zrb3 ; 

zmF : -pCzm , zrF3 ; 

zd: -sCsCzrCras/aa-rs/2 .03 , zl C los/oa-ls/2 . 03 3 , zcCaa»cas3 3 ; 
zo: -sCsCzr Cras/aa+rs/2 .03 , zlC las/aa*ls/2 .033, zcCaa*cas3 3 ; 
zs:-sCzrCrs3,zlCls33; 

Csolve For the network variables} 

denom: -csCcmCcaCzo, zmF , cl 3 , caCzo, zmb, cl 3 , cl 3 , cmCzd , zd, cl 3 , cl 3 ; 
vF : — cmplxC vs/2 .0, -vs/C2 . 0*a3 , cl 3 ; 
vb: -cmplxCvs/2 .0, vs/C2 . 0“a3 , cl 3 ; 

isF : -cdCcaCcmC vF , caCzo, zmb ,cl3,cl3,cmCvb,zd,cl3,cl3, denom, cl 3 ; 
isb: -cdCcaCcmC vb, ca Cza, zmf , cl3,cl3,cmCvF,zd,cl3,cl3, denom, cl 3 ; 

im:-caCisF, isb, cl 3 ; 
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ia : -cmCcmplxCO .0,3..0/a,cl3,csCisf,isb,cl3,cl3; 
i_line: “caCim, iq, cl 3 ; 

irf : -cdCcmCisf , zm , cl 1 , caCzm , zrf , cl 3 , cl 3 ; 
irb: -cdCcmC isb , zm, cl 3 , caCzm, zrb, cl 3 , cl 3 ; 

sef in: -cmCvf , ccC isf , 13 ,cl3; 
sebin: -cmCvb, ccCisb, 13 , Cl3 ; 
pern: -2 ,0»CrCsef in3+rCsebin3 3 ; 

mi_line: -sqrtCrC i_line3“rC i_line3*rC i_line*l 3 “rC i.l ine+13 3 
power .factor : -pein/ C vs*mi_line3 ; 

mirf : -sqrtCrC irf 3*rC irf 3+rC irf *13*rC irf *13 3 ; 
mirb: -sqrtCrC irb3*rC irb3*rC irb*13*rC irb+13 3 ; 
pf : -2 . 0*mirf*mirf *rr*C 1 . 0-slip3 /slip ; 
pb: -2 . 0*mirb*mirb*rr*C 1 .0-slip3/C2 .0-slip3 ; 
pm: -Cpf-pb3/pbase ; 

Cplot results! 

x: -5*raundCpm»sx3 ; 

y : - 190-round Cfm*sy 3 ; platCx,y,13; 

y : * 190-round C mi. line'sy /SO . 03 ; platCx, y , 1 3 ; 

y : “190-raundC power. f actor*sy 3 ; plotCx, y , 1 3 ; 

y : -190-roundCpm»sy*pbase/pein3 ; plotCx, y , 13 ; 

end; 

while not keypressed x:-l; 
end . 
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